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卷 积 神 经 网 络 方法 在 岛 礁 类 海啸 流水 动力 
特性 演变 的 应 用 


摘要 : 快速 预测 和 评估 海啸 是 预防 海洋 灾害 工作 中 的 重要 部 分 ， 对 于 海洋 工程 和 人 民生 命 财 产 安全 具有 
要 意义 。 本 文 以 一 维 卷 积 神经 网 络 (1-Dimensional Convolutional Neural Network, CON U 构建 岛 礁 
地 形 的 类 海啸 波 水 动力 特性 演变 模型 。 通 过 输入 类 海啸 波 波 高 时 程 曲 线 的 观测 值 ， 得 到 岛 礁 指定 地 点 的 水 
位 淹没 时 程 曲线 ， 实 现时 间 序 列 到 时 间 序 列 的 预测 ， 达 到 预防 海洋 灾害 的 目的 。 结 果 显 示 ， 对 类 海啸 波 到 
达 时 间 预 测 的 平均 误差 为 0.71%， 最 大 水 位 高 度 预测 的 平均 误差 为 6.99%，CONV1D 得 到 的 岛 礁 地 形 类 海 
啸 波 水 动力 特性 与 数值 结果 吻合 较 好 。 

关键 词 : 深度 学 习 卷 积 神经 网 络 海啸 预测 水 动力 特性 时 间 序 列 
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Application of Convolutional Neural Network Methods in the Evolution of 


Hydrodynamic Characteristics of Tsunamis like-wave over fringing reef 


Abstract: Rapid and accurate tsunami is an important part of Marine disaster prevention work, which is of great 
significance to Marine engineering and people's life and property safety. In this paper, based on 1-Dimensional 
Convolutional Neural Network (CONV1D), the evolution model of tsunami-like hydrodynamic characteristics of 
reef topography is constructed. By inputting observed values of wave heights resembling tsunami waves, the water 
inundation time series curves for specified locations on islands and reefs are generated. This achieves a prediction 
from one time series to another, serving the purpose of marine disaster prevention. The results indicate that the 
average error in predicting the arrival time of tsunami-like waves is 0.71%, and the average error in predicting 
maximum water levels is 6.99%. The hydrodynamic characteristics of island and reef terrains resembling tsunami 
waves obtained through CONV1ID exhibit a strong alignment with numerical results. 

Keywords: Deep learning, Convolutional neural network, Tsunami prediction, Hydrodynamic characteristics, 
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1 前 言 

中 国 南海 ,是 我 国 “建设 海洋 强国 ”战略 目标 的 重要 一 环 ， i ug LE 
轻重 的 意义 。 自 2013 FX, REO TA HET ASIA DUM, EE. E. Jm 
等 地 区 填 礁 造 陆 、 填 海造 陆 等 工程 活动 日 益 增 多 〈 李 庆功 等 ，2014; 姚 宇 等 ， iiis ， 由 于 该 
海域 位 于 亚 欧 、 太 平 洋 和 印度 -澳大利亚 板块 的 交汇 地 带 ， 包 含 马 尼 拉 海沟 、 苏 禄 海沟 和 苏 拉 威 
西海 沟 等 潜在 海啸 源 ， 因 此 风暴 潮 、 海 啸 等 海洋 极端 天 气 引 起 的 极端 波浪 事件 ， 是 工程 建设 必须 
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考虑 重要 因素 。 提 前 对 海啸 等 


p 
施 〈 包 澄 澜 ,2005) 。 


中 国 海 区 域 海啸 预警 中 心 主要 使 用 
传统 的 海啸 预警 方法 有 其 技术 方法 的 局 限 己 
造成 了 巨大 经 济 损失 和 人 员 伤 亡 〈 刘 烨 等 ， 
段 引起 人 们 重视 , 神经 网 络 在 此 之 前 已 初步 具 


2015; 


到 实现 。 然 而 传统 的 神经 网 在 海啸 预警 模拟 时 因 其 
实时 监测 的 海啸 预警 仍 是 
计算 机 算 力 的 快速 提升 ， 为 深度 学 习 发 


算 成 本 变 得 高 昂 ， 


然 灾害 进行 预测 和 评估 , 是 保护 人 民 的 4 


德国 SeisComP3 MEE 


EKRIR CER F, 


Antelope 地 震 监 ; 


徐 志 国 等 ， 


2022) 。 


E 命 和 财产 安全 重要 措 


测 处 理 系 统 ， 


E, 2011 年 日 本 东北 海啸 因 未 能 正确 评估 海啸 灾害 ， 
更 快 更 准确 的 海啸 预警 手 
备 预警 能 力 , 在 此 之 后 更 多 更 完善 的 神经 网 络 也 得 


庞大 的 计算 数据 和 复杂 的 非 线性 j 


d 程 , 使 得 计 


2021; Han et al, 


展 提 供 


了 必要 条 件 ， 


2021) 。 


2012 年 Hinton 小 组 发 明 的 深度 


卷 积 神经 网 络 (CNN ) AlexNet (Krizhevshky et al, 2012) 首先 在 图 像 处 理 方 面 取 得 成 功 ， 随 后 


深度 学 习 在 其 他 各 领域 引起 广泛 的 研究 。 
短期 记忆 神经 网 络 (LSTM) 得 到 广泛 应 月 
职 神 经 网 络 (CONVID) 在 很 短 的 时 间 内 对 心 


(Kiranyaz et al, 2015) 。 
2021 年 Makinoshima 等 


(Makinoshima et al, 2021) ， 它 区 别 于 传统 的 
2023)， 


时 刻 数 据 DE F, 


在 处 理 时 间 序 列 问题 上 ， 
H (Graves, 2012) 。 


循环 神经 网 络 


CRNNO 和 长 


2015 年 Kiranyaz 等 设计 出 一 维 卷 


应 用 程序 


有 图 进行 分 类 ， 成 为 了 最 先进 信号 处 理 
CNN-LSTM 混合 模型 可 以 解决 非 连续 数据 的 问题 〈 陆 继 翔 等 , 2019) 


构建 的 一 维 卷 积 


啸 发 生 时 ， 可 以 立即 、 


本 文 构建 时 间 序 列 至 
EEN AE i E 
列 变化 。CON1VD 由 于 


上 时 间 序 列 的 一 


地 形 的 传播 演变 数据 ,经 过 训练 精准 预测 岛 礁 礁 
其 较 低 的 计算 复杂 度 ， 


洋 工程 的 研究 提供 新 思路 科 


和 海洋 科学 的 交叉 发 展 ， 


2 模型 验证 


关于 海啸 数值 模拟 方法 ， 前 人 提供 了 许多 可 靠 上 
NHWAVE, HH Qu + (Qu et al, 2017; Qu et al, 2021) 
EIRIK, PREMK ERTE SEE 


促进 科技 创新 和 发 展 。 


只 神经 网 络 模型 使 
时 间 
实现 观测 数据 到 海啸 洪水 预 识 
准确 地 预测 指定 地 点 的 海啸 淹没 时 间 
维 卷 积 神经 网 络 海啸 预警 模型 , 使 用 数值 模拟 方法 得 到 类 
坪 的 某 一 位 置 海啸 来 临时 的 时 间 序 
适合 实 为 海洋 环境 和 海 


[新 方法 ， 拓 展 深 度 机 器 学 习 技术 的 应 用 


也 形 上 的 水 动力 学 过 程 。 


序列 。 


时 和 低 成 本 的 预测 ， 同 时 
领域 和 方 .推动 


提出 的 参数 化 类 海啸 波形 ， 


远海 观测 数据 进行 陆地 预警 
预测 方法 一 通过 前 几 个 时 刻 的 数据 推测 下 一 个 
1 的 时 间 序 列 到 时 间 序 列 预测 ， 


当 海 


机 器 学 习 技术 


方法， 本 文 是 基于 非 静 压 单 相 流 模型 


建立 高 精度 数 


参考 刘 等 (Liu et al,2021) 进 行 数值 的 


模拟 验证 。 实 验 数据 源 于 Roeber 等 (Roeber et al 2012) 在 海啸 盆地 建立 的 波浪 水 槽 物理 实验 ,我 们 


布置 的 二 维 
1 所 示 ， 左 侧 计 算 域 入 口 


水 槽 数值 模型 的 训 


1:5, WIEKEN 28 m, 


分 布 11 ANKAN 
小 dx-0.05 m; z 方 向 为 水 深 方向 ， 
水 深 为 1 m， 礁 坪 水 深 为 0m， 远 海水 深 为 1 m， 共 14 个 测 点 Roeber 
验证 结果 如 图 2 所 示 。 通 过 对 比 可 看 上 


为 造 波 边界 ， 右 侧 出 口 六 


| 算 域 布置 与 Roeber a 的 模型 布置 大 体 相 同 ， 


i m 
共 设 置 20 层 网 格 。 使 


上， 模拟 结果 与 实验 结 


x 方向 为 海中 波 传播 方向 


具体 布置 如 图 
前 斜坡 坡度 为 
网 格 大 


有 较 好 的 吻合 度 。 


j 非 静 压 模型 验证 入 射 波 高 为 0.5 m, 
等 的 不 透水 岛 礁 物理 实验 。 
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图 2 不 同 测 点 水 面 高 程 的 时 间 序列 (#8 =1.0m) 


Fig.2 Time series of water elevation at different wave gauges (H= 1.0m) 


3 数据 集 


图 3 为 基于 CONVID 的 类 海啸 波 神经 网 络 模型 的 示意 图 ,首先 使 用 NHWAVE 进行 了 117 


次 类 海啸 波 在 岛 礁 地 形 传播 演变 的 数值 模拟 ， 作 为 CONVID 的 数据 集 ,其 中 包括 106 组 训练 集 和 
11 组 测试 集 。 训 练 集 数据 是 模型 学 习 的 样本 ， 使 模型 能 够 从 输入 特征 中 学 习 到 数据 的 模式 、 关 


系 和 规律 ,从 而 能 够 在 未 见 过 的 数据 上 进行 泛 化 。 测 试 集 是 用 于 评估 机 器 学 习 模 型 性 能 和 泛 化 能 
力 的 一 组 数据 样本 的 集合 , 用 于 评估 模型 在 未 见 过 的 数据 上 的 表现 。 在 数值 模拟 中 ， 设 置 不 同 的 
初始 类 海啸 波 波 高 ， 最 小 波 高 0.20m 和 最 大 波 高 0.78m。 等 间隔 选取 初始 波 高 为 0.25m、0.30m、 
0.35m, 0.40m, 0.45m, 0.50m, 0.55m, 0.60m, 0.65m, 0.70m 和 0.75m 的 类 海啸 波 作为 测试 集 ， 


保证 其 具有 独立 性 、 代 表 性 、 标 签 或 目标 输出 、 足 够 数量 等 性 质 。 


测 指定 岛 礁 位 置 (WG11) 的 类 海啸 波 的 水 动力 特性 演变 过 程 ， 并 将 类 海啸 波 观 涡 
6s. 7s. 8s 和 9s 五 种 。 


E 


A 
Fr 
WE 


N p ex-sm sad 


观测 值 卷 积 层 平面 化 层 


图 3 卷 积 神经 网 络 结构 示意 图 
Fig.3 Convolutional Neural Network Architecture Diagram 


4 网 络 配置 


类 海啸 波 的 采样 频率 为 0.02Hz， 共 收集 了 16 秒 类 海啸 波 在 岛 礁 地 形 的 水 动力 特性 传播 演 
变 过 程 ， 使 用 10 个 位 置 CWGOI-WGI100 的 水 位 测 点 作为 CONVID 的 多 时 间 序 列 数据 输入 ， 预 


上 时 长 分 为 5s、 


CONVID 主要 由 十 二 层 神经 网 络 组 成 , 其 中 卷 积 层 包括 六 个 一 维 卷 积 层 和 三 个 池 化 层 组 成 ， 
最 后 是 三 个 全 连接 层 。 每 个 一 维 卷 积 层 之 后 都 有 一 个 Leaky Relu 激活 函数 层 (Maas et al, 2013), 
参数 Negative slope 为 0.01. 在 最 后 一 层 全 连接 层 前 加 入 随机 丢弃 层 ， 丢 弃 率 为 0.5， 表示 每 次 训 


练 输出 前 使 用 50% 特 征 数 进行 全 连接 ， 来 防止 过 度 拟 合 发 生 (Nitish et al,2014) o 


在 卷 积 层 中 ， 


设置 了 三 种 卷 积 核 ， 分 别 为 7、5、3， 移 动 步 长 为 1， 通过 填充 的 方法 使 得 卷 积 后 的 序列 长 度 不 
发 生变 化 ， 每 三 层 卷 积 层 后 设置 一 层 最 大 池 化 层 进 行 特 征 值 数量 的 优化 。 对 于 Ss, 6s. 7s. 8s 
和 9s 五 种 观测 时 长 ， 只 有 卷 积 后 得 到 特征 的 数量 有 所 不 同 ， 较 长 的 观测 周期 会 得 更 多 的 学 习 参 


数 见 表 1。 在 整个 训练 过 程 中 使 用 模拟 波形 和 预测 波形 的 均 方 误差 (MSE) 作为 损失 函数 ， 使 用 


Adam 优化 器 (Kingma et al，2014) 对 神经 网 络 中 的 权重 进行 优化 ， 将 损失 降低 至 


I 最小， 学 习 率 


为 0.0002, 其 他 参数 使 用 默认 值 .在 Pycharm 平台 上 基于 pytorch 框架 进行 模型 搭建 , 共 进行 8000 


次 训练 ， 使 用 i9-129000hx 处 理 器 和 4060 系 显 卡 完成 训练 大 概 需 要 60 分 钟 。 
5 训练 结果 


CONVID 训练 结束 后 ,我 们 使 用 测试 集 的 11 个 数值 模拟 结果 对 该 模型 进行 评测 。 表 1 为 不 


同 观 测 时 长 评估 类 海啸 波 的 模型 参数 量 和 需要 的 时 间 。 其 中 了 为 不 同 观 测 时 长 ,N 


F 


为 构建 模型 时 


116 ”所 需要 的 参数 量 ,Tr 为 预测 需要 的 时 间 。 由 表 1 可 知 对 于 较 长 的 观测 时 间 ， 需 要 使 用 更 多 的 参数 。 
117 Is 观测 值 的 增加 使 参数 的 数量 增加 了 25600, 主要 是 由 于 卷 积 层 后 全 连接 层 的 尺寸 较 大 ,其 训练 
118 时间 会 相应 增加 ,预测 所 用 时 间 也 有 增加 的 趋势 ， 实 际 时 间 由 计算 机 算 力 决定 。 使 用 CON1VD 进 
119 ， 行 海啸 评估 的 计算 时 间 只 有 不 到 0.1 秒 ,这 种 方法 比 传统 海啸 预警 机 制 和 数值 模拟 的 方法 的 快 得 
120 多， 而 且 需 要 计算 资源 更 少 。 从 上 述 数 据 可 知 ，CONVI1D 生成 岛 礁 测 点 水 位 时 程 曲线 的 效率 很 
121 高 ， 可 以 迅速 的 得 到 类 海啸 波 水 动力 特性 演变 过 程 ， 进 而 评估 和 预测 海啸 。 在 实现 中 ， 这 种 方法 
122 ”可 以 实现 对 海啸 等 极端 波浪 的 实时 监测 。 


=l 


123 
124 表 1 不 同 观测 时 长 模型 的 参数 量 和 预测 所 用 时 间 
125 Table.1 Number of parameters and prediction time used for different observation duration models 
观测 时 长 /Ts) 参数 量 /N 预测 所 用 时 间 /T。 © 

5 161110 0.0164 

6 186710 0.0175 

7 212310 0.0182 

8 237910 0.0201 

9 263510 0.0175 
126 
127 图 4 为 五 种 不 同 观测 时 长 时 程 曲线 结果 对 比 ， 其 中 图 4 Ca) 和 图 4 (b) 的 初始 类 海啸 波 波 
128 — 高 分 别 为 H-0.30m 和 及 =0.70m。 由 图 可 知 CONVID 预测 在 五 种 不 同 观测 时 长 下 的 结果 和 数值 
129 ”模拟 结果 都 可 以 较 好 的 吻合 ，=0.30m 的 类 海啸 波 随 着 观测 时 长 的 增加 ，WG11 测 点 的 水 位 的 
130 ”时 程 曲线 会 更 加 吻合 ， 有 =0.70m 的 类 海啸 波 随时 观测 时 长 的 增加 ， 出 现 到 达 时 间 提 前 和 最 大 水 


131 位 预测 不 足 的 问题 , 可 能 是 由 于 参数 的 增加 ， 导 致 模型 出 现 了 过 拟 合 问题 。 图 5 为 不 同 的 观测 时 
132 ”长 产生 的 误差 平方 和 (SSE) ， 由 图 4 和 图 5 可 知 ， 不 同 观测 时 长 模型 的 误差 平方 和 相差 较 小 ， 
133 ”在 5 秒 的 观测 时 长 下 可 以 达到 评估 类 海啸 波 水 动力 特性 的 要 求 ， 同 时 表现 出 CONVID 较 强 的 稳 
134 ”定性 和 泛 化 能 力 。 
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图 5 不 同 观测 时 长 的 误差 平方 和 


Fig.5 Sum of Squared Errors for Different Observation Durations 


图 6 HERA ISTE PAIS Fe] BE e E] PRU A EROR EE Eh E RTT 随 着 类 海啸 波 初始 波 高 变 大 ， 


测 点 水 位 的 时 程 | 


| 线 特征 特征 更 明显 、 特 征 值 更 多 ， 因 此 CONVID 计算 准确 度 更 高 ， 可 以 更 好 


的 展现 类 海啸 波 在 到 铝 礁 地 形 的 水 动力 特性 传播 演变 过 程 。 
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妈 6 六 秒 观测 时 长 下 不 同 波 高 预测 结果 对 比 


Fig.6 Comparison of Different Wave Height Predictions under a 6-second Observation Duration" 


为 了 进一步 评估 模型 预测 性 能 , 分 析 以 下 两 个 指标 : 类 海啸 波 到 达 礁 志 
的 最 大 水 位 。 选 取 五 种 观测 时 长 中 误差 


为 类 海啸 波 初始 波 高 、 


位 误差 、t op 为 wgll 观测 时 间 、t 


由 表 2 可 知 ， 该 模型 


n op ”为 WG11 观测 水 位 、n 


F 方 和 最 大 的 T-6s 模型 进行 分 析 ， 


p 


LE 


测 点 的 时 间 和 该 测 点 


如 表 2 所 示 ， 


为 WG11 预测 水 位 、E 


n 


为 预测 时 间 、 五 ， 为 时 间 误 差 。 


KH 


为 水 


现 的 最 高 水 位 的 最 大 误差 为 14.46%、 最 小 误差 为 2.64%、 最 高 水 位 


平均 误差 为 6.99%; 到 达 时 间 最 大 误差 为 1.709%、 最 小 误差 为 0%， 平 均 误差 0.71%。 图 7 为 不 


同 初始 波 高 的 最 大 水 位 和 到 达 时 间 误 差 散 点 图 , 由 图 
动 较 小 且 误差 很 小 ,两 者 都 呈现 出 有 规 得 


可 知 最 大 水 位 的 误差 波动 较 大 , 时 间 


误差 波 


的 波动 。 通 过 表 2 和 图 7 可 知 ， 类 海啸 波 的 水 动力 特 特 


信 度 很 高 。 


性 难以 精准 的 预测 ， 需 要 构建 复杂 度 更 高 的 CONV1D 模型 ， 但 


是 该 模型 


K 2 六 秒 观测 时 长 下 不 同 波 高 的 最 高 水 位 和 到 达 时 间 


Table.2 Maximum Water Level and Arrival Time of Different Wave Heights under a 6-second Observation 


Hím f] g/m N, /m 
0.25 0.0550 0.0527 
0.30 0.0616 0.0592 
0.35 0.0706 0.0629 
0.40 0.0878 0.0751 
0.45 0.0843 0.7066 
0.50 0.0872 0.0836 
0.55 0.0897 0.0849 
0.60 0.1018 0.0884 
0.65 0.1046 0.1016 
0.70 0.1049 0.0991 
0.75 0.1025 0.0964 
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Fig.7 Error in Maximum Water Level and Arrival Time for Different Wave Heights under a 6-second Observation 


Duration 
6 结论 


本 文通 过 构建 CONCID 模型 ， 实 现 类 海啸 波 在 岛 礁 地 形 水 动力 特性 演变 过 程 的 预测 。 该 模 
型 得 到 的 类 海啸 波 预 测 结果 和 数值 模拟 结果 基本 吻合 ， 对 于 不 同 观 测 时 长 得 到 的 结果 基本 一 致 。 
在 最 大 的 误差 平方 和 的 情况 下 ， 模 型 评估 岛 礁 指定 地 点 的 最 大 水 位 误差 平均 值 为 6.99%， 到 达 时 
间 误 差 平 均值 仅 为 0.71%。CONV1D 只 需要 较 少 的 计算 资源 获得 较 高 的 计算 效率 ， 可 以 提前 精 
准 获取 波浪 信息 ， 达 到 预防 灾害 和 保护 人 民生 命 财产 安 的 目的 。 
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